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formula for event horizons. 
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I. INTRODUCTION 



Fundamentally, there are two ways to characterize a black hole. The first focuses on 
causal structure and says that a point in an asymptotically flat spacetime is inside a black 
hole if no signal from that point can reach future null infinity (see, for example, [1]). The 
boundary of the black hole region is the event horizon. This is an intuitive definition but 
at the same time is teleological and so highly non-local — one must trace all causal curves 
from a point before deciding whether or not it lies inside a black hole. By this definition 
neither large spacetime curvatures nor singularities are necessary for black hole existence. 

In contrast, the second characterization is quasilocal and geometric, saying that a point in 
spacetime is inside a black hole if it lies on a trapped surface. Such surfaces, from which both 
future directed null expansions are everywhere negative, are indicative of large spacetime 
curvature and so this definition directly makes the association between strong gravitational 
fields and black holes. It is also directly connected to spacetime singularities as the existence 
of a trapped surface is sufficient to imply the existence of a spacetime singularity [2].^ 

Traditionally, the second viewpoint inspired the definition of an apparent horizon. Given 
a Cauchy surface, the trapped region is defined as the (closure of the) union of all the 
trapped surfaces contained in that slice of spacetime. The apparent horizon is then the (two- 
dimensional) boundary of the trapped region [1] . Correspondingly, if a region of spacetime 
is foliated by Cauchy surfaces, then one can locate the apparent horizon on each slice and 
so define a time-evolved (three-dimensional) version of the apparent horizon. Often this is 
also referred to as the apparent horizon. 

For any given foliation of a spacetime, most trapped surfaces will not lie in the specified 
slices. Thus the time-evolved apparent horizon is defined by only a subset of the total num- 
ber of trapped surfaces and so is certainly slicing dependent and contained in the "total" 
trapped region. The time-evolved apparent horizons defined by various foliations will typ- 
ically intersect each other multiple times and also will usually have fully trapped surfaces 
lying partially outside of them (see for example [3-5] for discussions on these points). 

By definition, an apparent horizon is a boundary between regions containing trapped and 
untrapped surfaces. As such, it is no surprise that it is marginally outer trapped — that 
is the expansion, of its outward null normal, i, vanishes [1]. Now, while it is certainly 
not true that all such surfaces will be apparent horizons, in practical terms it is clearly 
easier to find the marginally outer trapped surfaces (MOTS) and then identify apparent 
horizon candidates, rather than try to proceed by first identifying all trapped surfaces. This 
is the approach taken in numerical relativity (see for example [6]) and in fact in that field 
"apparent horizon" is usually understood to mean the outermost MOTS. 

Other properties can also be expected from an apparent horizon. First if there are 
fully trapped surfaces "just inside" the apparent horizon, then there should exist arbitrarily 
small inward deformations generated by some spacelike vector field X"- under which the 
outward expansion becomes negative, i.e. SxG{e) < 0. Furthermore, by continuity with the 
fully trapped surfaces, the expansion 6'(„) of the inward-pointing null normal n should be 
negative on the apparent horizon. Over a decade ago, Hayward formalized this intuition in 



Of course the two definitions are not entirely independent. In spacetimes with a well-defined future null 
infinity, trapped surfaces necessarily lie within the causally defined black hole region. Further the two 
characterizations both identify the same region for the family of Kerr-Newmann solutions [1] (though this 
is not necessarily true in more general spacetimes). 
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his definition of trapping horizons [7]. Here, we will be mainly interested in future outer 
trapping horizons (FOTHs) which are three- surfaces that may be foliated by two-surfaces 
with 9(^£) = 0, 6'(„) < and 6n0(i) < 0. Although these will include most time-evolved 
apparent horizons, the focus has now shifted to the three-dimensional horizon surface itself. 
There is no reference to a foliation of the full spacetime. 

Marginally trapped surfaces are practical as they may be (relatively) easily identified 
in simulations as well as studied with standard geometric tools. However, they are also 
philosophically appealing as, in contrast to event horizons, their evolution is causal. As such, 
this idea of studying the "boundary" of a trapped region without first finding a corresponding 
bulk has become increasingly popular in the last few years. Apart from the further studies 
of trapping horizons [8, 9] there are closely related programmes such as isolated horizons 
which identify and study equilibrium states [10-14], dynamical horizons which correspond 
to non-trivially evolving horizons [15-25], and slowly evolving horizons which are "almost 
isolated" trapping horizons [26, 27]. For reviews of the field see, for example [28-30]. 

Trapping and dynamical horizons each come equipped with a preferred foliation into two- 
surfaces. It is this foliation which is used to define the null normals and hence the expansions. 
In this paper we will focus on variations of these two-surfaces as a route to a better under- 
standing of the existence, evolutions and deformations of the full three-dimensional horizons. 
Given a vector field X = Ai — Bn (where A and B are functions) normal to a horizon cross- 
section, the corresponding variations are generated by 5x- Particularly important is the 
variation SxO{£) which turns out to be a second-order elliptic operator in B that is defined 
by the intrinsic and extrinsic geometry of the two-surface along with components of the 
Einstein tensor. The techniques used are similar to those described in [18, 31], however the 
emphasis is somewhat different. In those papers, the focus was on two-dimensional MOTS 
in a three-dimensional slice of space-time, whereas we consider general deformations of the 
two-surface in the full four-dimensional spacetime. 

Solutions of 5x0(11) = generate both the evolution and the possible deformations of 
FOTHs. A judicious application of a maximum principle to the resulting elliptic partial 
differential equation is a key to both deriving new results about these horizons and also 
unifying much of the existing knowledge under a common formalism. Among other results, 
this technique will be used to show that any foliation of an isolated FOTH may be freely 
deformed (a well-known result) while the foliation of dynamical FOTH is rigid (a related 
version of this was first shown in [16]). Conversely we will see that an isolated FOTH is rigid 
against normal deformations (it may only be deformed into itself) while this is definitely not 
true for a dynamical FOTH. However, the allowed deformations in the dynamical case are 
strongly restricted by rules that are consistent with, although slightly different from, those 
seen in [16]. 

Results from black hole physics also follow from the deformation equations. Apart from 
the second law of horizon dynamics [7, 15, 18] and angular momentum flux laws [15, 23, 32] 
we will examine slowly evolving horizons in some detail [26, 27]. In particular we will explain 
how a horizon may be invariantly characterized as "almost" isolated. From this definition we 
will examine the circumstances under which a FOTH has a well-defined and slowly varying 
surface gravity and derive the first law for slowly evolving horizons. The related flux laws 
for event and dynamical horizons also follow from the variation equations. 

The plan of the paper is as follows. We begin, in Section II, by considering the geometry 
of general spacelike two-surfaces embedded in four- dimensional spacetimes and study how 
that geometry changes if the surfaces are deformed . From there we specialize to two-surfaces 
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that satisfy 6'(£) = 0, < 0, and < and study the properties of deformations that 

preserve those conditions. This is done in Section III. Next, in Section IV, we apply these 
results to gain a better understanding of isolated, dynamical, and future outer trapping hori- 
zons. With this general understanding in hand we turn to a more specialized study of slowly 
evolving horizons in Section V. Finally we compare the flux laws for slowly evolving hori- 
zons to the corresponding ones for dynamical and event horizons in Section VI. Numerous 
technical results are found in Appendices A-C. 



II. TWO-SURFACES AND THEIR DEFORMATIONS 



To begin, we review the differential geometry of two-surfaces embedded in four- 
dimensional spacetimes. Most of the results appearing in this section are not new but 
it is useful to gather them together here both for reference and to set the notation and 
emphasis that will be found in future sections. 



A. Two-surface geometry 



Let S* be a closed and orientable two-surface that is (smoothly) embedded in a four- 
dimensional time-oriented spacetime (M, gab) which has metric compatible covariant deriva- 
tive Va. Then there are just two future-pointing null directions normal to S. Let and 
n"- be null vector fields pointing in these directions; in situations where this is meaningful 
we will always take £° and as respectively outward and inward pointing. If we further 
require that i ■ n = —1 then there is only one remaining degree of rescaling freedom in the 
definition of these vector fields. 

The intrinsic geometry of S is defined by the induced metric and area-form. The definition 
of these quantities is independent of the choice of null vectors above, however for our purposes 
it is most useful to express them in terms of these vectors. Thus, the induced metric on S 
can be written as 



Qab = gab + Lnb + hria , (2.1) 

while the area two- form e satisfies e = € A n A e where e is the four- volume form on M. This 
metric also defines the compatible covariant derivative operator da and (two-dimensional) 
Ricci scalar R on this two-surface. 

The extrinsic geometry describes how S is embedded in M and in the usual way is defined 
by how the (in this case null) normal vectors change over S. The extrinsic curvatures are: 

k'-2 = q:qi^M and = g^g^V^n, . (2.2) 



These are symmetric since (.a and are, by 
way we decompose them as 

^ab = \dmab + (Tab and 

where the expansions 

Oil) = r'^aib and 



definition, surface forming. In the standard 
A^l? = l0^n)qab + al? . (2.3) 

0in) = r'^anb, (2.4) 
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are the traces of the extrinsic curvatures and the shears 

<y^I = {qlqi - \qa,t^ V da and ai? ^ {qlq^ - \qa,t^ V.n, , (2.5) 

are the trace-free parts. 

The last part of the extrinsic geometry is given by the connection on the normal cotangent 
bundle T^S namely 

cu, :=-g^n,V,r. (2.6) 

To see that this is the connection consider a general normal one-form = aia + (3na G T^S. 
Then a direct calculation shows that 

g^Ve/if, = /c^^^ + (4« + ^aaYb + (4/5 - ^al3)nb , (2.7) 

where k^J^^ = ak^^^ + (3kl^\ Thus, the covariant derivative on this normal bundle is 

{aib + (3nb) := (4^ + Cjaa)ib + (daP - 0JaP)nb , (2.8) 

and Ua clearly acts as the connection. The gauge dependence in this case is the scaling chosen 
for the null vectors. If £ — > and n ^ n/f for some function /, then the corresponding 
transformation for the connection is 

Ua-^'^a + dalnf . (2.9) 

As usual the geometric, gauge independent, information associated with the connection is 
contained in its curvature which in this case is 

Qab = daUJb - dbUa , (2.10) 

and this is constrained by the four-space curvature via the Ricci equation (Appendix A 1): 

^ab = q:<ittnfC,.ef + a^M:^ - Vi^) , (2.11) 

where Cede/ is the Weyl curvature of the full spacetime. In this paper we will usually be 
more interested in the connection itself rather than this curvature. 

Other constraints relating the geometry of S to the full four-space curvature come from 
the Gauss and Codazzi equations. The Gauss equation is 

€<il&nef,H = Rabcd + {k^IktJ + fci^^S) - (fcl?e; + kt'Jkf,) , (2.12) 

where Tlefgh and Rabcd are the (four- and two-dimensional) Riemann tensors, while the 
(slightly modified) Codazzi equations are 

{da - Ua)9it) = 2(4 - Ub^J^ - q'aGbct - itCbcdeti^n^ and (2.13) 
(4 + ^a)Bi^n) = 2(4 + Cob)alr^'' - tGbcu' + 2g^acden'=£V , (2.14) 

where Gab = T^ab — ^T^dab is the Einstein tensor. A derivation of these relations can be found 
in Appendix A. 
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curve of constant s 

FIG. 1: A schematic of a section of 7$ around So- The variation vector field X° is everywhere 
tangent to the tube, points along curves of constant s, and generates the foliation. Thus for 
sufficiently small e, one can intuitively write = So + eX. 

B. Deforming a two-surface 

1. Defining variations 

A variation or deformation of a two-surface So is a smooth, one-to-one function A) : 
So X [— Ao, Ao] — > M (with Aq some real number) such that $(S'o, 0) = So- Thus, $ generates 
a (finite) three-surface 7$ and that surface is foliated by images Sx = $(S'o,A) of So as 
depicted in Fig. 1. The variation vector field X"" = (d/dX)"" is tangent to the curves of 
constant s G So- The flow generated by this vector field maps leaves of constant A into each 
other. Unless otherwise noted, we will restrict our attention to normal variations where X" 
is everywhere perpendicular to the Sx and so can be written: 

X'' = Ar-Bn\ (2.15) 

for some functions A and B- There are no restrictions on the values of A and B- However, 
in later sections we will usually assume that X°- is i-oriented so that A > 0- Then if S > 0, 
X°- is spacelike while B < means that it is timelike. We will mainly be interested in 
situations where B > 0, hence the negative sign in (2.15). Fig. 2 should help to keep the 
various cases straight. 

For any values of A and B, the map $ deforms So into successive surfaces Sx- To quantify 
the change in the geometry, we identify points along the curves of constant s and then 
calculate derivatives of the intrinsic and extrinsic geometry with respect to the parameter 
A. Taking the normal connection Ua as an example, we write its variation as Sx^a and note 
that calculating this quantity amounts to: 

1. using $ to construct the Sx in a neighbourhood of A = 0, 

2. constructing the Ua on those Sx (among other things this will involve choosing a scaling 
of the null vectors), and 
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Future 




FIG. 2: The normal space at a point s G 5a- Normal vectors written as X°' = Al"" — Bn"" point 
into the causal future if ^ > and B < Q. The shaded region represents the ^-oriented vector 
fields for which ^4 > and B may take any value. 

3. calculating the Lie derivative C x^a and puUing-back the result onto So- 
Thus, once the two-surfaces are constructed and uja calculated, we have 

5xuJa = ql^l^xi^c, (2.16) 
in standard four-dimensional notation. 

2. Calculating variations 

We now calculate some of these variations for the geometric quantities that are of interest 
in this paper. The easiest calculation is the variation of the two-metric. A couple of lines of 
algebra shows 

5x~qa, = Ak^2 - = (Ae^e) - 5^(n))ga. + 2(AalJ - EaS^) , (2.17) 

from which it follows that 

5xe = iAe^e)-B9^r.))e. (2.18) 

These expressions justify referring to 6'(£) and 6'(„) as expansions and o"^^^ and a^^^ as shears; 
respectively they describe how So expands and shears if deformed in the null directions. 

Finding the variations of the extrinsic quantities is more involved. Here we just outline 
the calculations but more details can be found in Appendix A 2. First we note that since 
£a and are everywhere normal to the Sx, both = and 6xna = (with the usual 
pull-backs understood). Thus 

X'VbL = -{da-Coa)B + Kxia and (2.19) 

X^VkUa = {da + Cda)A- Kxna (2.20) 
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where 



(2.21) 



is the component of the connection on the S\ normal bundles in the X-direction. Explicitly, 
under rescalings £ — > and n ^ n/ f , 



Kx K.X + Cx'^'^f ■ 



(2.22) 



We will usually refer to kx as the surface gravity associated with X"- in analogy with the 
corresponding quantities on a Killing or isolated horizon (though at this stage we make no 
claims about the physical content of this nomenclature). 

The importance of (2.19) and (2.20) is that they allow us to convert expressions involv- 
ing derivatives off of So into ones containing only quantities defined on So plus the gauge 
dependent surface gravity k^. Then, with the help of these relations a direct calculation 
shows that 



Sx(^{e) 



KxO. 



daCo" - RI2 + Gabtn^ 



-A [llaWf + G,,r£^+(l/2)^2j 



(2.23) 



where ||a)|p = cDacD" and ||cr*^^)p = af^a^^^"-^. 

We will also need to know SxO{n) and this is most easily calculated by substitutions into 
the above expression (2.23). Exchanging ia and Ua and sending A —B and B —A it 
is straightforward to see that 



dau'^ - R/2 + Gabny" - 6^)6 



5x0{n) + ^^xO^n) = d'^A + 2Cj''daA + A 

+ i? [11^^"^ f + Gafen V + (1/2)^2^J . 

Finally, the variation of the normal connection (also derived in Appendix A 2) is 



6xC0a - daKx = 'k^^ [d' A + A] - k'^^ [d^^ B - Cj' B] 



(n) 



(2.24) 



(2.25) 



where r*^ = Ai'^ + Bn'^ is normal to X'^. Enlisting the help of the Codazzi equations by 
combining A x (2.13) + B x (2.14) — 2 x (2.25) this can be rewritten as 



dx^a + iAOf^i) - BO^n 



da^x - d, {Aa^^ , + Ba^^^ J + q^G^cT'' 

+ -^da {A9(^i) + B9(^n)) — G{i)daA — 9(^n)daB , 



(2.26) 



which eliminates the Weyl dependence. This is the form that we will use. 

These variations will be sufficient for most of our considerations. Note that equivalent 
or closely related versions of the expressions for 6x0{e), 5x6{n), and 5x^a have previously 
appeared in, for example, [3, 18, 23, 25, 30, 33], though not with the particular two-surface 
emphasis that we adopt here. 
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C. Angular momentum and its evolution 



Physically, the connection Ua defines the angular momentum associated with any rotation 
vector field 0" on a closed two-surface S [9, 15, 19, 23, 26, 32]. By definition the fiow 
associated with such a 0" has only two fixed points and foliates the remainder of S into 
closed integral curves of parameter length 27r. The canonical example of such a field is a 
Killing vector field of the two-metric Qab and in that particular case it replaces the fiat-space 
notion of an axis of rotation. However even if it is not a Killing vector field, it is standard 
to define the angular momentum of S relative to 0" as 



= ^ / er^a . (2.27) 

s 

Note that any rotation vector field is necessarily divergence free and a quick calculation 
with the help of (2.9) shows that this expression is independent of the scaling of the null 
vector fields. Alternatively if the surface S has a suitable topology such as S*^, then a 
divergence free 0" can necessarily be written in the form 0" = e°'^dh( for the area-form iab 
and some function (. Then, the scaling independence is made explicit if we rewrite (2.27) 
as 

where as usual we drop the indices and write in bold any form that is being integrated over. 

We now consider how the angular momentum changes under deformations. To this end 
we multiply (2.26) by the area form e and so obtain 

Sx {eCua) = e [daKx - dh (Aaf ' + Ba^:^') + q^'Ghcr') (2.29) 

Now, extend off of S by demanding that 6x4'"' = — essentially this is equivalent to 
the flat space requirement that angular momentum be measured relative to a flxed axis of 
rotation. Then, contracting (2.29) with 0" and integrating over S, several total divergence 
terms vanish and we flnd that 

SxJ[<P] = ^j^^ {Gabrr' + a(^):a(<^) + A</.«4^(,) + , (2.30) 

where the shear with respect to a general vector fleld w"" is 

^if = {qlqi - Iqam'") "^cw, , (2.31) 

and a^'^'^-.a^'^^ = a^^Ja^'^^"^ (this double contraction notation is adopted from [23]). 

We will return to equation (2.30) in SectionVC where we will be able to neglect the 
last two-terms and so interpret the change in angular momentum as coming from a flux 
of stress-energy (with the help of the Einstein equation) and a flux of shear. A detailed 
discussion and fluid mechanical interpretation of (2.29) and (2.30) can also be found in [23]. 
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D. The "constraint" law 

Finally, before specializing to the two-surfaces associated with horizons we derive the 
following relation. First, combining A x (2.23) + B x (2.24) we find that 

KxOx = A6x0^e) + B6x0^n) + da{A<rB - Bd^A - 2AB0'') (2.32) 



2" WW ' 

where we again have = Aia + Bua and 6(^x) and 6'(t-) are defined in the obvious way. 
Integrating this over S, the total divergence term vanishes and we find that 



(2.33) 



More generally for a vector field 

JY'^ = + x'' , (2.34) 
where x°' is everywhere transverse to the 5* we can combine (2.29) and (2.32) to obtain 

^ / {nxSxe + i^Sxii^^a)} = ^ f e [G.fcA'V^ + a^^^ia^^^] (2.35) 
sttG Js ovrG Js 



1 



SttG 
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A6x9ie) + B 5x6 in) + ^^(x)^(r) 



In cases where kx is constant and x° = Vtcf)"" for some constant Vt (which is not related 
to the curvature of the normal bundle) and rotation vector field satisfying 5x(t>'^ = 0, the 
left-hand side of this equation takes a particularly familiar form: 

'*^-d + fij[0]. (2.36) 



The similarity to the first law of black hole mechanics is not coincidence. In sections V 
and VI we will see that the dynamical version of the first law for both event and trapping 
horizon are closely related to (2.35). 

A version of this relation was referred to as the horizon constraint law in [17, 26] due to 
its equivalence to the integrated diffeomorphism constraint on 7$. A discussion of its terms 
and their interpretation if the horizon is viewed as a viscous fluid can also be found in [23] . 
Furthermore, in [24], equation (2.35) is interpreted as a second order evolution equation for 
the area. 



III. FUTURE OUTER TRAPPED SURFACES 

The expressions of the previous section hold for any spacelike two-surface embedded in 
any four-dimensional spacetime. In this paper however, our main interest will in the two- 
surfaces that foliate future outer trapping horizons which in turn are embedded in solutions 
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of the Einstein equations. Thus in this section we consider spacehke two-surfaces on which 
Gab = 8iiTab and for which 6'(£) = 0, 9(^n) < and there is a scahng of the null vectors such 
that 6nd(e) < 0. Adapting Hayward's nomenclature we will call such surfaces future outer 
trapping surfaces (FOTS). 

First, we consider the conditions under which a marginally trapped surface S with 6'(^) = 
and d^n) < is a FOTS. To this end, we we set A = 0, i? = —1 and = in (2.23) and 
so find that on a marginally trapped surface 

5J^^^ = -R/2 + \\uf - dau" + {STxG)Tabtn^ ■ (3.1) 

Then it is immediate that the 5n0[t) < condition is determined entirely by this component of 
the stress-energy tensor along with the intrinsic and extrinsic geometry of - no derivatives 
need to be taken off of the surface. That said, checking this condition is slightly more 
complicated than just calculating this quantity with an arbitrary scaling of the null vectors. 
For example, in Appendix C it is shown that for the standard scaling of null vectors on a 
Kerr horizon, 6n0(e) is not always less then zero. A rescaling is necessary for this relationship 
to become apparent. 

Now, for any spacelike two-surface on which = 0, equation (2.23) simplifies to become 
5x0ii) = -d^B + 2u''daB - B5Jf^,) + A5,e^,) , (3.2) 
where 5n0i^t) takes the form shown in equation (3.1) and 

= -ll^^'^f - {^TiG)Tabtt . (3.3) 

If this is a FOTS we can adopt a scaling so that 5n0[t) < 0, while if we assume the null 
energy condition it also follows that 6e6(^i) < 0. Then, we can draw several conclusions about 
FOTS. 

First strengthening to the dominant energy condition, there is the well-known [34] re- 
striction on the topology of such a surface. 

FOTS Property 1 If S is a closed and orientable FOTS on which the dominant energy 
condition holds, then it is homeomorphic to S"^. 

This follows from equation (3.1). On integrating this over 5* and doing a bit of rearranging 
we find that : 

X = ^ [ H-Snen) + Ipf + STrTabtn'} , (3.4) 
^vr Js 

where x is the Euler characteristic of S. Now, by assumption 6n0{£) < while the dominant 
energy condition implies that the matter term is positive. Thus x > and this is sufficient 
to tell us that 5* must be homeomorphic to a two-sphere since that is the only closed and 
orientable two-surface with positive Euler characteristic. □ 

Next, we consider the conditions under which a FOTS may be deformed whilst preserving 
its defining characteristics. Sufficiently small variations will always leave 6'(„) < and 
^n^ie) < and so the key to understanding these deformations is finding normal vector fields 
X" such that SxO{e) = 0. We assume that all fields are at least twice differentiable. 

We begin with the case where = everywhere on a FOTS. Then we have: 
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FOTS Property 2 If S is a FOTS on which S£6(^£-^ = everywhere, then variation vectors 
X"- satisfy 6x0(e) = and only if they are parallel to i"". 

Starting with the first part, if = everywhere on S then equation (3.2) becomes 

5x9^e) = d'^B - 2u''daB + BSJ^^^ , (3.5) 

and it is trivial that B = ^ SxO{e) = for any value of A. It is also straightforward to see 
that the converse must be true. First applying the maximum principle of Appendix B to 
^xO{e) = we find that B must be either constant or everywhere negative. Similarly applying 
the corresponding minimum principle to —5xd{t) = we find that B must be either constant 
or everywhere positive. Thus, B must be constant and it is clear that since 7^ that 

constant must be zero. The result is established. □ 

It is also true that all such deformations leave the intrinsic geometry of S invariant: 

FOTS Property 3 Let S he a FOTS on which 5^9 i^i) = everywhere and the null energy 
condition holds. Then 6x0{e) = implies that SxQab = 0. That is, all deformations leave the 
intrinsic geometry invariant. 

If the null energy condition holds then all terms in (3.3) are non-negative and so if (5f6'(£) = 

they must all, including cr^^\ be zero. Then with i? = by Property 2 and 6'(£) = by 
assumption, equation (2.17) implies that SxQab = as required. □ 

Such results are familiar from the isolated horizon literature and we will return to them 
in section IV. For now however we consider FOTS on which Se6(£) is somewhere non-zero. 
In doing this we restrict our attention to ^-oriented variation vector fields for which A > 
(Fig. 2). 

FOTS Property 4 Let S be a FOTS and assume the null energy condition. Then, if 
^e^ii) 7^ anywhere on S, all i-oriented variation vectors X°' that satisfy SxO(i) = are 
spacelike everywhere on S . 

If the null energy condition holds then 6£6(^i) < by equation (3.3). Thus with A > 0, 
equation (3.2) implies that 

-d^B + 2cu"45 - BSn9(i) > , (3.6) 

everywhere on 5". Then by the minimum principle of Appendix B, B is either everywhere 
positive or everywhere constant. If it is constant then the derivatives in (3.2) vanish and 

everywhere on 5*. In particular this must hold at the point where < and so with 

A > and 5„6'(£) < we again find that B > 0. Thus in either case gabX"'X^ = 2AB > 
and X°- is spacelike. □ 

Combining this with Property 2 we see that FOTS satisfying the null energy condition 
may be cleanly split into two classes - those for which ^-oriented variation vector fields 
satisfying SxO{e) = are null and those for which these vectors are spacelike. Such an X" 
cannot be timelike and what is more it cannot be partly null and partly spacelike. 

We also know something about how the geometry of a FOTS must change with respect 
to a spacelike 9(^£) = preserving deformation: 
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FIG. 3: Given two spacelike variation vector fields X" and X'°' on a FOTS that satisfy 6x0{e) = 
5x'd(e) = 0, then if they are not constant multiples of each other they must interweave so that X'°' 
will sometimes point into the future of X"' and sometimes will point into the past. 

FOTS Property 5 Let S be a FOTS and assume the null energy condition. Then if 5(,9(^t) 
is non-zero anywhere on S, 6x^ > everywhere. The deformation causes S to expand 
everywhere. 

If 7^ anywhere on S, then by Property 4, X"- must be everywhere spacelike with 

B > 0. Then 6x€. = —B9(^n)e where —B6(n) > since 6'(„) < by assumption. □ 

Finally, it is quite clear that in general there will be an infinite number of X°- that will 
solve 5x0(t) = and so an equally infinite number of FOTS-preserving deformations. For 
example, if 5i9{^t) is nowhere zero and we choose any B G C^{S) then we can always solve 
5x0{i) = to find a corresponding A (though unless B is constant, there is no guarantee that 
the resulting X" will be ^-oriented). Furthermore if 5^6' = everywhere, then Property 2 
tells us that for any X" = At, 5x0{i>) = 0. 

In contrast to the = case where all allowed variation vector fields must be parallel 
to t, if 5^6* (£) is somewhere non-zero then (apart from constant rescalings) no two variation 
vector fields are parallel. Instead they must interweave as shown in Fig. 3. 

FOTS Property 6 Let S be a FOTS with 6i6(^i) somewhere non-zero and assume the null 
energy condition. Further let X*^ and X'"- be two i-oriented, FOTS-preserving deformation 
vector fields. Then either 

1. X'" = AoX" for some constant Ao or 

2. X'" interweaves X" in the sense that X'°'Ta takes both positive and negative values on 
S , where r'^ = At + Bn"' is the usual forward pointing timelike normal to S. 

By Property 4, both X" and X"^ are spacelike and so with X" = At — Brf" for some positive 
A and -B, we have 

X'" = (a + \)At - XBrf (3.8) 

for some functions a and A that satisfy both A > and a + A > 0. 

Now, to begin our analysis, let us consider the case where a takes both positive and 
negative values. Then we immediately see that 

X'Va = -aAB , (3.9) 
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and so this corresponds to the second posited behaviour for X''^. 

By contrast, if a does not take both positive and negative values then at least one of 
a < or a > must be true. As preparation to exploring these two possibilities, we note 
that given 5x0(i) = 0, it is straightforward to see that Sx'0(£) = reduces to 

B(fX + 2(rf"5 - BCj'')da\ = aAdeOf^i) . (3.10) 

Thus, if a < this equation implies that 

B(fX + 2{d''B-BCu'')da\>0. (3.11) 

Now, A must achieve a maximum on S, so we can define a new function A* = A — Xmax 
that also satisfies equation (3.11) and achieves a maximum value of zero. Then, by our 
usual maximum principle, since A* is not everywhere negative, it must be constant. This 
means that A is also constant and so by (3.10), « = 0. Thus a < ^ X'" = XqX"' for some 
constant Aq. Similar reasoning shows that a > also implies that X'" is a constant multiple 

IV. HORIZONS AND THEIR PROPERTIES 

In this section we apply the properties of FOTS to gain a better understanding of future 
outer trapping horizons. First though we recall some definitions. 

A. Horizons 

We begin with the definition of a future outer trapping horizon (FOTH). A trapping 
horizon is a three-dimensional submanifold of a spacetime (M, Qab) that may be foliated with 
closed and spacelike two-surfaces Sy (where f is a foliation parameter) on which ^(^) = 
[7]. Trapping horizons are classified by the values taken by 6'(„) and SnO(e) on their leaves. 
A trapping horizon is said to be future (past) if 6'(„) < (6'(„) > 0) while it is outer (inner) 
if there is a scaling of the null vectors such that (5„,^(£) < {6n6{e) > 0).^ These names are 
taken from the horizons that satisfy these conditions in a fully extended Kerr or Reissner- 
Nordstrom spacetime. Thus the event horizon is a future outer trapping horizon (FOTH) 
and the inner Cauchy horizon is future inner. The corresponding white hole horizons are 
past outer and past inner respectively. In this paper, we will be mainly in interested in 
FOTHs (which are fohated by FOTS). 

As noted in the introduction, apart from (future outer) trapping horizons, there are other 
closely related quasilocal horizons. Here we recall the definition of isolated and dynamical 
horizons, while in section V we will consider slowly evolving horizons. 

A three-dimensional submanifold of a spacetime {M,gab) is a non-expanding horizon if: 
i) it is null and topologically S* x IR for some closed two-manifold S, ii) 6(^e) = and iii) 



^ The original definition of [7] is phrased in terms of a dual-null foliation of the spacetime in some vicinity of 
the horizon instead of the variations that we use here. In constructing such a foliation, one must usually 
abandon the i ■ n = —1 normalization [33] , however having done this the definition can be phrased in 
terms of Lie derivatives rather than variation operators. That said, the definitions are equivalent. 
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—T°'^ii, is future directed and causal [11]. As usual is an outward pointing normal and we 
note that since the horizon is null, no foliation is required for its construction. However, a 
foliation is certainly no hindrance to a three-surface being a non-expanding horizon, and any 
null FOTH satisfying the null energy condition will certainly be a non-expanding horizon. 

Non-expanding horizons are the simplest objects in the isolated horizon family. Any 
non-expanding horizon can be turned into a weakly isolated horizon if the scaling of the null 
vectors is chosen so that C = for 

UJa := -ribV a^t = -KiUa + C^a , (4.1) 

where the arrow indicates a pull-back into the cotangent bundle of the non-expanding hori- 
zon. With this scaling zeroth and first laws of isolated horizon mechanics may be established 
[11]. Furthermore, an isolated horizon is obtained by strengthening the above conditions 
to require that the entire extrinsic geometry encoded in the derivative operator be time 
independent. 

Finally, a three-dimensional sub-manifold of a spacetime (M, gab) is a dynamical horizon 
if it: i) is spacelike and ii) can be foliated by spacelike two-surfaces such that the null normals 
to those surfaces satisfy 9(^t) = and < [15]. As we shall see in the next subsection, if 
the null energy condition holds and 6i9{^t) 7^ then dynamical horizons are FOTHs and vice 
versa. 

Non-expanding, isolated and dynamical horizons each have a rich set of properties that 
may be derived directly from their definitions. However, examples of spacetimes exist that 
contain (and are actually foliated by) isolated horizons but do not contain trapped surfaces 
[14]. Similarly there are spacetimes with dynamical horizons but no trapped surfaces [35]. 
Thus if we take trapped surfaces as the defining property of black holes, neither of these 
definitions is sufficient to specifically single out black holes. Instead they represent necessary 
conditions which must be supplemented to become sufficient. 

Given this observation, we will take FOTHs as our basic objects and classify them by 
hybridizing the naming systems. Thus a FOTH that also satisfies one of these other sets of 
properties will be referred to as non-expanding, weakly isolated, or dynamical as appropriate. 

B. Properties of FOTHs 

A FOTH can be thought of as the variation surface associated with a (finite) deformation 
of a FOTS (like the one shown in Fig. 1). Specifically given a FOTH and a foliation labelling 
f , we can always find a tangent vector field V that is normal to the and which satisfies 

Cvv = l. (4.2) 

Then, V" can be viewed as a variation vector field and as for other deformation vectors we 
write 

for some functions A and B. Further, identifying points on the different by the flow 
generated by V, we write derivatives of the two-geometry with respect to v as 

l^S.. (4.4) 
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Then we may apply our results on FOTS to learn about FOTHs. 

First, their topology is strongly constrained [7] and this follows directly from FOTS 
Property 1: 

FOTH Property 1 Let H be a FOTH and assume the dominant energy condition. Then 
H has topology S"^ x IR. 

Next, we consider the circumstances under which a FOTH is non- expanding and those 
by which it is dynamical. In particular we will be interested in transitions between these 
behaviours and so as a preliminary we define non-expanding and dynamical sections of a 
FOTH. If a FOTH H is null (that is B = everywhere) for some range f i < f < ^2 of the 
foliation parameter then we will refer to this as a non-expanding section of H. In contrast 
if H is spacelike (that is AB > everywhere) for some range v[ < v < v'2 of the foliation 
parameter then we will refer to this as a dynamical section of H. 

Then by FOTS Properties 2 and 4, on any of a FOTH either B = everywhere or 
B ^ anywhere. Thus, no element of the foliation is partly non-expanding and partly 
dynamical. Transitions between non-expanding and dynamical sections must happen "all- 
at-once" . 

FOTH Property 2 Let H be a FOTH with foliation and assume the null energy con- 
dition. Then H may be completely partitioned into non- expanding and dynamical sections. 
On non-expanding sections 5i9(^i) = everywhere while on dynamical sections 6e6(^i) < at 
least somewhere on each S^. 

This surprising result was first shown (using slightly different assumptions) in [18]. As 
one application, it guarantees that a FOTH on which 6i9(^£) 7^ for at least one point on each 
cross section is necessarily a dynamical horizon. The converse is slightly more restricted and 
we must require that 6e6(^£) 7^ everywhere to ensure that a dynamical horizon also be a 
FOTH. To see this consider that a dynamical horizon is necessarily spacelike and so X"' can 
always be oriented so that A > and -B > 0. Scaling the null vectors so that B is constant 
on each cross-section of the horizon (3.2) simplifies to 

Bdne^e) = A5ee^£) . (4.5) 

Thus, if S£9(^i) 7^ everywhere on a cross-section then the null energy condition guarantees 
that is strictly negative and so 5„,^(f) < everywhere as well. Consequently such a dynamical 
horizon is a FOTH. However, if 6e6(^i-) = somewhere (or everywhere) on the surface, then at 
those points 6n0(i) will also vanish and the dynamical horizon will not be a FOTH. Examples 
of such dynamical horizons are considered in [35] . 

Next by FOTS Properties 3 and 5, dynamical versus non-expanding sections differ in 
more than just signature: 

FOTH Property 3 Let AH be a (section of a) FOTH and assume that the null energy 
condition holds. 

1. If AH is non- expanding, then the intrinsic geometry of the is invariant. That is 
SvQab = 0. 

2. If AH is dynamical, then Sy locally increases in area everywhere. That is 6ye = Fe 
for some function F that is everywhere positive. 
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The first part of this property is well known from the isolated horizon literature [11, 12]. 
The second part is essentially Hayward's second law of black hole mechanics [7]. FOTHs 
expand if and only if they are dynamical. Otherwise their intrinsic geometry is unchanging. 

Having reaffirmed these basic results, we can consider how FOTHs may be deformed 
while preserving their defining characteristics. First we consider variations of the foliation 
that leave H itself invariant. 

FOTH Property 4 Let H he a FOTH and assume that the null energy condition holds. 

1. If H is non-expanding then the foliation may he smoothly deformed hy any vector field 
of the form X" = where f is any positive function. 

2. If H is dynamical then the foliation is rigid and can he relahelled hut not deformed. 

The first part of this property follows directly from FOTS Property 2 which we apply to 
all of the simultaneously^. The second part follows from FOTS Property 6 which tells 
us that if (5v^(^) = then 5{^aV)d{i) = if and only if a is a constant over Sy. For such a 
constant a the variation would preserve (but relabel) the foliation. 

Thus, the foliation may be deformed in an infinite number of ways on a non-expanding 
(section of a) FOTH while it may not be deformed at all on a dynamical section. The first 
part of this property is consistent with the fact that particular foliations are not important 
for most of the isolated horizon formalism. The second part is consistent with the more 
general result of [16] which says that if 5* is a FOTS and is a dynamical FOTH with 
foliation S^, then S (Z H ii and only if it 5* = for some v. 

Next, we consider the other ways in which a FOTH may be deformed. Again the non- 
expanding and dynamical cases are quite different. 

FOTH Property 5 Let H he a FOTH and assume the null energy condition. We consider 
variations that deform H hut preserve the FOTH conditions. 

1. If H is non-expanding then all allowed variations map H into itself. 

2. If H is dynamical then all allowed variations deform each leaf of its foliation partly 
into the causal future and partly into the causal past of H . 

This property also follows directly from FOTS Properties 2 and 6 and can be rephrased 
to say that non-expanding FOTHs are stable against deformations but dynamical FOTHs 
are not stable in his way. The second part is consistent with [16] where it was shown that if 
Hi and H2 are dynamical FOTHs with foliations S^^ and 6*^2 respectively, then no S^^ can 
lie entirely in the causal past (or future) of H2 and vice versa. Thus, either S^^ is causally 
disconnected from H2 or it intersects it so that it lies partly in the causal future and partly 
in the causal past. Property 5 may be thought of as a local version of that global result. 



^ Such variations of the foHation will generally only be allowed for a finite range of the variation. Beyond 
that range we may violate one of the defining conditions, for though the intrinsic geometry of H will remain 
invariant, the extrinsic geometry of the two-surfaces normal to H will change. Thus, finite variations may 
ultimately generate foliation slices with or Sn9[i) non-negative. 
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V. SLOWLY EVOLVING HORIZONS AND THEIR FLUX LAWS 



A significant part of pliysics deals witli systems tliat are at or near equilibrium. For 
liorizons, we naturally take non-expanding FOTHs to be equilibrium states. Included in 
this class are the event horizons of all (non-extreme) Kerr-Newmann black holes (whose 
properties are summarized in Appendix C). Further, such FOTHs are automatically isolated 
horizons and so all of the results from that formalism apply to them. Recall too that by 
FOTH Property 5, these surfaces (locally) can only be deformed into themselves, so there 
is no ambiguity about their exact location. 

It is then natural to consider near-equilibrium horizons, which should be those that change 
slowly in time. Now, while this is a very reasonable condition intuitively, it is not so easy 
to geometrically and invariantly characterize such horizons — keep in mind that dynamical 
FOTHs are spacelike and so there is no natural notion of time intrinsically associated with 
the surface. In this section we will motivate a definition of these slowly evolving horizons 
and then explore some of its consequences. A version of the definition and many of the 
properties appeared in [26] but most of the derivations appear here for the first time. 

A. Slowly expanding horizons 

Intuitively, we would expect the properties of a slowly expanding horizon to be "close" 
to those of a non-expanding or isolated horizon. In this section we will formalise this 
requirement. This is essentially done by performing a perturbation expansion around a non- 
expanding horizon in powers of a small parameter e. However, as we shall see, care must 
be taken as spacelike and null surfaces are fundamentally different, and in particular it is 
non-trivial to find a normalization which transitions smoothly between the dynamical and 
isolated regimes. 

In order to proceed, we will slightly modify the formalism introduced over the last three 
sections. First we loosen it so that while the evolution vector field V" still generates the 
foliation, it is no longer tied to a specific labelling; that is we only require that C\;v = a{v) 
for some positive function a{v). Next we tighten it by choosing A = 1 so that 

V^ = r-Cn% (5.1) 

where we change notation from i? to C to avoid confusion with the more general variations 
considered previously. The net effect is to reduce the scaling freedom of the null vectors to 
that contained in a{v) so that for a fixed foliation labelling 

= aV" , r = , = -< and C = c?Co , (5.2) 

a 

where the subscript "o" indicates a quantity defined with a = 1. 

If a FOTH is null and non-expanding, then C = and the evolution vector field V = I. 
Thus, the most direct way to define a slowly expanding horizon might seem to be to say that 
it is a (section of) a horizon on which C is very small. Unfortunately this is not a viable 
strategy as one can use the rescaling freedom (5.2) to make C arbitrarily small on any AiJ. 
Instead we proceed by focusing on how the are evolved by V". 

For our purposes non-expanding horizons have two key properties: they are null and the 
intrinsic geometry of their foliating two-surfaces is invariant. Our invariant characterization 
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of a slowly expanding horizon draws on both of these ideas. First, we consider the evolution 
of the area-form e on S^. From (2.18) we have 

Cye = -C9(^ri)^, (5.3) 

which is certainly scaling dependent. However this dependence may be easily isolated by 
rewriting 




£vi=-C^(n)6 = ||V|| -W-^(„)e , (5.4) 



so that the rescaling freedom is restricted to | |V| |. The term in parentheses then provides an 
invariant measure of the rate of expansion. Among other properties it vanishes if the horizon 
is non-expanding and on a dynamical horizon section is equal to the rate expansion of e with 
respect to the unit-normalized version of the evolution vector field. We will consider the 
rate of expansion to be slow if this is small and this approach is borne out by the examples 
of [27]. In particular if —\fC9{n) is of order e and the scaling of the null vectors is chosen so 
that ||V|| is commensurate (as would be reasonable for an "almost-null" surface) then 
will be of order e^. 

This notion of a slow area change captures the essence of a slowly expanding horizon 
however we still need to require both that the surface be "almost-null" and that the rest 
of the intrinsic geometry also be slowly changing. We have already seen that restrictions 
on the norm of are not invariant, so instead we will implement these ideas together by 
requiring that the evolution of the two-metric be characterized by the expansion and shear 
associated with (as they would be for a truly null surface). That is from (2.17), 



Further, the matter flow across the horizon similarly should be, to lowest order, the same 
as that across a null surface. That is we expect 

V'TS^Tabti", (5.5) 

where Ta = ia + Cua is the usual timelike normal. Given that these quantities vanish on 
a non-expanding horizon we would also expect both of these quantities to be small. An 
invariant implementation of these ideas gives rise to the following definition: 

Definition: Let AH be a section of a future outer trapping horizon foliated by spacelike 
two-surfaces Sy so that AH = {UyS^ '■ vi < v < V2}- Further let V" be an evolution vector 
field that generates the foliation so that C^v = a{v) for some positive function a(w), and 
scale the null vectors so that V" = — Cn°'. Then AH is a slowly expanding horizon if the 
dominant energy condition holds and 

1. e <^ 1 where 

e^/R]j = Maximum [C (Ha^^^f + TabU^n'' + ^(„)/2)] , (5.6) 

2. and Tabi°"n^ ^ l/Rj^, where Rh is the areal radius ^ya/Air of the horizon, 
and 
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3. two-surface derivatives of horizon fields are at most of the same order in e as the 
(maximum of the) original fields. For example, ||(iaC|| ~ Cmax/Ru, where Cmax is 
largest absolute value attained by C on S^. 

Throughout this definition and what follows, an expression like X "^Y means that X < k^Y 
for some constant ko of order one, and in particular k^, ^\/e. 

As we will shortly see, the first condition is simply a scaling invariant way of writing our 
earlier requirements that the horizon geometry be slowly changing and that those changes 
be dominated by the ^-components of quantities. Additionally, it introduces a simple way 
to define a value for e on each cross section of the horizon. If the a{v) is chosen so that 

C ^ e^ (5.7) 

then this scaling of the null normals is compatible with the evolution parameter e. 

The second and third conditions are restrictions on the horizon geometry. Effectively, 
they ensure that the geometry of the horizon is not too extreme. These conditions are 
rather mild, and as a partial justification for these assumptions in Appendix C it is shown 
that they all hold on a Kerr horizon with the standard foliation. 

We now consider the implications of these assumptions. To begin, by equation (3.1) our 
restrictions on the horizon geometry immediately imply that 

\Sn9(e)\^^. (5.8) 

Note the use of the absolute value sign. Even though we must have Snd(e) < for some 
scaling, in general this need not be true for the particular A = 1 scaling that we have 
chosen. 

Next, we can bound the flux of matter and the gravitation shear at the horizon. To do 
this, we make use of the fact that (5v^(£) = 0. Then, equation (3.2) along with the magnitude 
of C fixed by (5.7) implies that 

\5ee^e)\ = Ik^'^f + 8nGT^,ti' < , (5.9) 



or, using the null energy condition 



ll^^'^f ~ -BT and MGT^.t^' ^ ^ ■ (5.10) 

Thus these terms are bound by the size of e. Explicit examples of these bounds in action 
may be found in [27] where both matter and shear driven expansions are considered in some 
detail. 

From these results it follows that the two-metric is slowly changing with the highest order 
contributions coming from quantities associated with We have 

C vqab = 2a^^ + (Ce^n)qab - 2Cal;:^) , (5.11) 
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while 



Cvl= -C9^n)e . (5.12) 
o(.2) 

In addition, there is a bound on the rate of change of the intrinsic curvature. The evolution 
of R is derived in (A20) and is given by 

CvR = 2d^d\a^^ - Ca^^) + d\Ce^n)) + Ce^n)R • (5.13) 

From this we see immediately that C\)R ^ ^/R?h- Thus, the intrinsic geometry of the 
horizon is slowly varying, at a rate e. The choice (5.7) effectively scales V"" to reflect the 
slowly expanding nature of the horizon. In a transition to isolation so that e — > this choice 
will force ||V|| — > 0, thus ensuring that the limit is continuous. 

We can also bound the flux of energy momentum through the horizon. We have already 
seen that Tabi°'i'' ^ e^/Rf^. Furthermore, this is the main flux of energy through the horizon 
as our definition implies that 

SnGV'Tj^n = STiGTabti^ - SyrGC^I^^riV . (5. 14) 

0(e2) 0(e4) 

If the dominant energy condition holds, then there is a further constraint on the components 
of the stress-energy tensor. In that case, with r° future-directed and timelike on a dynamical 
FOTH, — TafoT* must also be future directed and causal. Then, 

g''\T^cr%ndT'')<0 \\qtTbcr'f<2iTatrT')iT,,n'^T''). (5.15) 

However, by (5.6), (5.7) and (5.10), {Tad^r^T^dn^r'^) is of order / R\j. Thus, 

llfcri < ^ . (5.16) 

This result, in conjunction with (2.13) can be used to bound one of the components of the 
Weyl tensor. We obtain 

g^acde«V<^. (5.17) 
Rh 

This is equivalent to \&i = Cabcd^°'T^'^^'^n'^ ~ ^/^"h- isolated horizon the equivalent 

quantities vanish, namely q^Tbci'^ = and \E'i = [11]. 

The flux of incoming gravitational radiation is encoded in \E'o = C abcd^""'^'^ ^'^^'^ ■, another 
of the components of the Weyl tensor. For a horizon to be slowly expanding, one would 
expect this quantity to be small. We can show that this is the case by considering the 
evolution of the shear a^^\ derived in (A15). Keeping only the lowest order terms, we obtain 

C vaiJ - Kvai? = &'&^C,def + 0{^) . (5.18) 

Therefore a'^^'^ will remain small only if \E^o ~ ^1 R^h- 

Finally, we turn to = Cabcd^'^TTT'^fn'^n'^- Begin by noting that \E'2 is invariant under 
rescalings of i and n. On an isolated horizon, the value of \&2 is not restricted, however, 
£v^2 = 0. For a vacuum, slowly evolving horizon, £v^2 ^ ^/^h- This follows from the 
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Gauss and Ricci relations of Appendix A 1 which can be used to rewrite \l/2 in terms of the 
intrinsic and extrinsic horizon geometry. The result then follows directly from the fact that 
< e/ Rh and the fact that cD, R and a^^^ are all slowly evolving. 
To summarize, we have seen that the definition of a slowly expanding horizon captures 
many expected properties of a near equilibrium black hole. Specifically, the intrinsic geome- 
try, including the area and two-curvature is slowly changing and there is little flux of matter 
or gravitational energy through the horizon. Of course, the given orders of quantities are 
bounds rather than requirements on the size of those terms. For example, on a spherically 
symmetric slowly expanding horizon, cr^^ vanishes identically and so the metric is unchang- 
ing at first order. Similarly in a vacuum spacetime where a horizon grows through the 
absorption of gravitational waves, all matter terms will vanish. Examples of both of these 
behaviours may be found in [27]. 



B. Slowly evolving horizons and the first law 

In the previous section, we have placed requirements on the intrinsic geometry of the 
horizon and arrived at the notion of a slowly expanding horizon. Now, we shall impose some 
further restrictions in order to obtain the first law of black hole mechanics. This will be done 
by restricting the extrinsic geometry of the horizon. The inspiration for the extra conditions 
comes from the weakly isolated horizons defined in section IV A. A non-expanding horizon 
becomes weakly isolated if the scaling of the null vectors is chosen so that 

C,Ua = Q. (5.19) 

In this case it is also true [12] that one can always find a "good" foliation of surfaces such 
that 

C,e^n) = Q. (5.20) 

With this foliation, a suitable scaling sets = — [(if]^ and (5.19) can be decomposed as 

CiKi = Q and £ ^u)a = . (5-21) 

We will enforce versions of these conditions to obtain slowly evolving horizons. However, 
there is a distinction between the slowly evolving and isolated cases. Since non-expanding 
horizons do not have a pre-determined, fixed foliation, it is always possible to rescale the null 
normal ^ so that conditions (5.19) and (5.20) are satisfied. In contrast, a non-equilibrium 
horizon comes endowed with a unique foliation, so the equivalent conditions are not guar- 
anteed to hold, they will have to be checked. This motivation leads us to the following: 

Definition: Let Aif be a slowly expanding section of a FOTH with a compatible scaling 
of the null normals. Then it is said to be a slowly evolving horizon if in addition 

1. ||£v'^a|| and |£v'^v| ~ ^/R^h 

2. \Cve (n) I ~ 
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The first consequence of tlie above definition is that on a slowly evolving horizon, the 
surface gravity is slowly varying. It follows immediately from the definition that is slowly 
changing in time.^The fact that it is nearly constant across each two-surface follows from 
(2.29). Keeping only the lowest order terms, we have: 

£ v^a = daKv - d'a^^ - MGq%,t + 0{^) . (5.22) 
From the definition above, (5.10) and (5.16), it follows immediately that 

Ma«^v|| ~ -4- • (5-23) 

That is, the surface gravity is approximately constant over each slice of the foliation. Since 
we have also required the surface gravity to be slowly changing up the horizon, it follows 
that over a foliation parameter range on A/J that is small relative to 1/e, 

= + 0(e) (5.24) 

for some constant Kq. Note however, that if a FOTH is slowly evolving for long enough, 
then larger changes can accumulate. 

Slowly evolving horizons obey a first law of black hole mechanics [26]. Applying the 
slowly evolving horizon conditions to equation (2.33), it reduces to a dynamical version of 
this first law. To order we obtain: 



SttG Js I SnG 

where a is the area of the two-surfaces and d = £ye. Interestingly, in [24], equation 
(2.33) has been interpreted as a second order evolution equation for the horizon area. In 
the slowly evolving limit, this reduces to (5.25). In section VI we will compare this version 
of the first law other well-known flux laws, but for now consider it in its own right. 

Let us examine the two energy flux terms contributing to the area increase. The first is 
the square of the gravitational shear at the horizon, while the second is the flux of matter 
stress energy through the horizon. The first term is necessarily positive, and provided the 
null energy condition holds, so is the second. Then, FOTH Property 3 implies that for a 
dynamical slowly evolving horizon Ko > - the average surface gravity of a slowly evolving 
horizon is necessarily positive. See [36] for further discussion of this point and its relation 
to extremal horizons. 

Next, we would like to examine whether (5.25) can be integrated to give a value for the 
horizon energy. For a slowly evolving horizon, (2.24) reduces at leading order to 

KoOi^n) = -R/2 + STTGTabtn^ + + daC^" + 0(e) . (5.26) 

For a spacetime which is close to spherically symmetric, such as those considered in [27], 
the Ua terms are of order e or smaller while if the only matter is radially infalling dust the 
matter term may also be neglected. Then 

R 1 



^ As for a weakly isolated horizon [11] this condition can equivalently be implemented by imposing conditions 
on the various quantities that arise if one takes the £v derivative of equation (2.24) and then deriving 
the desired result for Cynv- However for simplicity we just impose the condition directly. 
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If one scales the null vectors so that = — 2/i?//, which is the value taken in the 
Schwarzschild spacetime, it follows that Kq = 1/2Rh- Then it is immediate that (5.25) 
may be integrated to give an energy oi E k, Rh/2. In this case, the first law can be written 
as 

where we have taken the foliation label v to be compatible with null scaling (that is £ yf = 1). 
Thus we recover all of the standard notions of black hole mechanics: the energy is given 
by the Smarr formula and its time rate of change may be written in terms of both KqCI and 
a flux law. In more general situations however, things are not quite so tidy. While (5.25) 
always holds, away from spherical symmetry and in the presence of alternative matter fields 
the later simplifications cannot be made. Thus, in general it is not guaranteed that (5.25) 
will integrate to a tidy expression for the energy - this is not too surprising given the well- 
known uncertainties in defining (quasi)localized gravitational energy. 



C. Approximate symmetries and angular momentum 

We would like to generalize the first law for slowly evolving black holes to include angular 
momentum. We begin by noting that on a slowly evolving horizon, (2.30) simplifies to 

r [- (£). (0) -^ 

where both terms are 0(e) due to (5.10) and (5.16) respectively. Thus, the angular mo- 
mentum associated to any rotation vector field must be slowly varying, even if is not a 
symmetry of the horizon. However, in this case the change in the angular momentum is only 
restricted to be at most of order e while the area (and energy) evolve at a rate proportional 
to e^. This reflects the fact that we have not required the vector field to be a symmetry of 
the horizon. Since the horizon is not in equilibrium, we do not expect it to possess an exact 
symmetry, and instead introduce the notion of an approximate symmetry. 

Definition: Let /\H be a section of a FOTH and cf)^ G TS^, be a rotation vector field as 
defined in section II C. Then 0" is said to be an approximate symmetry of the horizon if: 

2. \C^d(.n,)\ < e/R]j, and 

3. |87rGr,50V^| < e^/Rl. 

The first two conditions require 0" to be an approximate symmetry of the intrinsic and 
(part of) the extrinsic geometry of the two-surfaces. The third condition is not really a 
symmetry condition but instead says that the angular matter flux should be particularly 
small in the 0" direction (in general (5.16) only restricts such fluxes to be of order e). These 
conditions are sufficient to guarantee that the angular momentum measured relative to 0" 
changes to order with the expression for the rate of change given as before by (5.29). 
Furthermore, the change in angular momentum is proportional to a gravitational plus a 
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matter flux. As for other slowly evolving relations, these fluxes are calculated as if /\H was 
a null surface. Finally, we note that the condition Cycj)"' = condition can be weakened to 
||>Cv0"|| ^ e^/R^fj, without affecting the angular momentum evolution. This allows for slight 
changes in the approximate symmetry direction as the horizon evolves. 

Now, let us consider the first law for slowly evolving, rotating horizons. In this case, we 
allow for a more general evolution vector C,"", tangent to the horizon but with components 
both normal and tangent to the cross sections. We restrict the allowed in the following 
manner. First C, should preserve the foliation of the horizon. That is ^'^ = + for some 
normalization of V" and some which is tangent to the horizon cross sections. Second, this 

should generate rotations. That is, it should integrate to a flow which foliates the cross 
sections into two fixed points plus a congruence of closed curves and further those closed 
curves should have a common period (see [17] for a further discussion of rotations). Finally, 
we require that respect the slowly evolving nature of the horizon in a non-trivial way. 
That is we require that ||,^|| ^ e while the norms of C^qab, '^^^(n) and C^uJa should be of 
order e. Then, we can write 



where 0" is an approximate symmetry and Vt is an angular velocity which is constant on 
each horizon cross section. ^ In addition, we require that VL changes only very slowly with 
respect to V so that ~ I^h- The evolution vector field ^ is then the analogue of 

the Killing evolution vector field on a Kerr horizon. Thus, combining (5.25) and (5.29) and 
expanding to order we obtain the first law for rotating horizons: 



To this order the fluxes can be calculated using ^ + Vtcj)"-. 

As in the last section this familiar form of the first law will hold for all horizons of the 
considered class. Again however, care must be taken in the choice of scaling for the null 
vectors and choice of Vt if one hopes to be able to integrate it to give a simple energy expres- 
sion for a rotating FOTH. That said, for perturbations around a Kerr horizon it is possible 
to scale the null vectors accordingly and so obtain the standard functional dependence of 
energy on area and angular momentum. 

VI. COMPARISON WITH OTHER FLUX LAWS 

The forms of the first law derived in the last section are not new. Other dynamical flux 
laws may be found in the literature which apply to other types of horizons. In this section 
we see how two of these laws may also be derived from the deformation equations of section 
II and compare them to our form of the first law. 

A. The first law for event horizons 

A version of the first law usually known as the Hawking-Hartle formula can also be derived 
for event horizons. We now outline how it arises with our chief focus being a comparison 



^ Note that the angular velocity Vl is not related to the curvature of the normal bundle ^ab- 



a 



(5.30) 




(5.31) 
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with the first law for slowly evolving horizons. For further details (and a slightly different 
derivation) see the original paper [37]. 

An event horizon is a null surface ruled by a congruence of geodesies. Thus its evolution 
is governed by Raychaudhuri's equation [1] which from our point of view is either equation 
(2.33) or (3.2) with A = 1 and B = 0: 

SiOi^t) = ^Ai) - h^'\? - Gabtt - (1/2)^?,) . (6.1) 

While this is not strictly necessary for a null surface, we'll assume that the horizon is foliated 
by spacelike two-surfaces and that the foliation is compatible with affine scalings of the 1°'. 
Thus, we could scale this null evolution vector so that = 0. However, even if we don't, it 
is immediate that is a function of v alone and so dai^^ = 0. 

Then, multiplying by the area element and writing variations in Lie derivative form, we 
have 

+ ^(,)/2) Cee = Ce{Cee) + e{ f + Gabti') ■ (6.2) 

If we think of this as a perturbation of a standard (stationary) black hole solution then the 
scaling of the null vectors should be such that the surface gravity is of order 1 / while 
other quantities appearing in the above equation should be close to zero. Thus we should 
have 9{^i) -C i^i and so can drop the 9(£)/2 on the left-hand side of the above equation. Then, 
integrating both over and "up" the horizon between two surfaces Si and 5*2 we find that 

/ dv{K^a) ^ d\ll + dv e{\\a^^^f + Gabte^) ■ (6.3) 

J Vl J VI J 

While it is certainly true that an event horizon always expands until it reaches an ultimate 
equilibrium state, it is also true that during periods of quiescence when nothing much is 
happening it can approach an isolated horizon. In particular, as discussed in [30], d can 
become arbitrarily close to zero. If we consider an evolution between two such "equilibrium" 
states, the first term on the right-hand side of (6.3) can also be neglected. Thus, we arrive 
at the Hawking-Hartle formula: 

l'V2 pV2 p 

dv{Kid)^ / dv / e(||a(^)f + G„,,r£'') . (6.4) 
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This is the event horizon analogue of our (5.25). However, the equation does not imply 
a causal relationship from fluxes to changes in area. This reflects the teleological nature of 
event horizons. Expansions of event horizons are caused by an absence of interactions, while 
fluxes through them instead force decreases of the rate of expansion (again see [30]). One 
of the significant implications of the Hawking-Hartle formula is that time averages smooth 
out these strange behaviours. 

Despite the very different character of FOTHs and event horizons, there are remarkable 
similarities between the first laws for slowly evolving horizons (5.25) and event horizons (6.4). 
In both cases we must impose a condition which forces the horizon to be slowly evolving in 
order to obtain the first law. Thus, they both only hold near equilibrium. Further, neither 
of these laws either specifies (or requires) an energy definition for the black hole. 

If we further assume that the event horizon has an approximate symmetry generated by a 
spacelike vector field 0" so that Cecp" and C^O(^i) can be neglected then (2.30) again becomes 



s 



/ e{^-^ + T^,r<p'\ . (6.5) 
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In contrast to (6.4) this is a snapshot rather than time- integrated flux law. This is a reflection 
of the much broader applicability of the angular momentum flux law which holds not just 
for horizons but for any surface [15, 17, 23, 32, 38, 39]. 

Then, given a slowly varying angular velocity Q, we can define an evolution vector field 

= r + n(f)'' , (6.6) 

and combine (6.4) and (6.5), to get the more general dynamical first law for event horizons 

Apart from the time-integration this is, of course, the same as the corresponding law (5.31) 
for slowly evolving horizons. 



B. Dynamical horizon flux law 

Finally, we consider the dynamical horizon flux law of [9, 15]. Despite the by now familiar 
flux terms that appear in this equation, it is different from the first laws that we have already 
seen. Instead of linking a na term to fluxes of gravitational and matter stress energy through 
the horizon, it is concerned with how these fluxes change the energy associated with the black 
hole. As we shall see, it is not clear how the approaches can be directly connected. 

Let if be a FOTH and further let E{v) be any function that is increasing whenever H is 
dynamical but remains constant whenever (if ever) it is isolated. The obvious choice is some 
function of the area but in principle other any other function with this property would work 
just as well. Thanks to FOTH Property 2 (which excludes the possibility of FOTS that are 
partially isolated and partially dynamical) we can always scale the null vectors so that 

B dE , ^ 

where G is the gravitational constant. This means that B will be constant on each leaf of 
the foliation®. With a bit of rearranging, equation (3.2) integrates over to become 



dE_ 

dv 



e 

Sv 



(6.9) 



where we have used the fact that the S^, is topologically S"^ to rewrite eR = Sn. Then, 
integrating up the horizon we find that 



ff2 

E{v2) - E{v,) = I dv{ I ~e 
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(6.10) 



If E{v) is any state function of the horizon (such as energy or entropy) this is interpreted 
as a flux law. This flux law is valid for any functional E{v), although Hayward [9] has ar- 
gued that the Hawking (or irreducible) mass is the most natural choice. Given the dominant 



This scaling is equivalent to requiring that the pull-back of £a to T*H satisfy J_ — —^dv. 
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energy condition, the right-hand side of (6.10) consists of three non- negative terms. The 
interpretation of the first two is reasonably straightforward as a matter flux and a flux of 
gravitational energy through the horizon. However, the third term is not so easily under- 
stood. A possible interpretation is that this represents a flux of rotational energy [9, 15, 19] 
however this seems unlikely as uJa is associated with the angular momentum itself rather 
than its flux. 

A direct comparison in the slowly evolving limit with (5.25) cannot be made due to 
the different methods of scaling the null vectors. The slowly evolving formalism can be 
generalized to allow for such a comparison by relaxing the requirement that A = 1, however 
even then the limit does not go through directly. To see this note that we can expand the 
integrand of the right-hand side of the dynamical horizon flux law as 

A(||aWf/(87rG)+T,,r£'') +5(||cu||V(87rG')+T,,£V) . (6.11) 

In order to make a comparison with the slowly evolving horizons, we consider the limit 
where A ^ 1 and B ^ e^. Then, all of these terms are of order and the last two — which 
do not arise in the slowly evolving law — cannot be neglected. Further, for the apparently 
common terms, the A ^ 1 version of the slowly evolving constraint law (2.33) gives rise to 
shear and flux terms 

(||(TWf/(87rG) + Tahiy') . (6.12) 

These differ by a factor of A from those in (6.11). 

Clearly in the limit where the horizon is slowly evolving, this flux law will reduce to (5.25) 
up to the factor of A discrepancy discussed above. Now, it is quite possible that as a horizon 
approaches equilibrium, A will be constant at leading order over each cross-section. If this is 
the case, the two results will agree as we can set A = 1 -|- 0(e). However, there does not seem 
to be any analytic justification for this, and the examples of slowly evolving horizons studied 
so far in [27] are all (approximately) spherically symmetric, whence A is automatically 
(approximately) constant by construction. Thus, whether or not this assumption holds is 
currently an open question. 

Motivated by the results for slowly evolving horizons, we can derive an alternative flux law 
which is valid on all dynamical horizons. On a dynamical horizon, we must have Se6(^£) < 0, 
and SnO(^i) < with strict inequality at some point on the horizon. If these conditions do not 
hold, then the horizon will not be spacelike. Making use of these conditions, we can rewrite 
(3.2) as 

-BSnee = -ASeO^i) . (6.13) 

On integration over a cross section of the horizon, both sides are guaranteed to be positive. 
Therefore, making use of (3.1) and (3.3) we obtain 

(1 - e) f ) = / ^TT^ + AW (6.14) 

where 



2G / ./c 87rG 



Cj'\\ + 2GTai,t'rf (6.15) 



and e is related to the extremality of the horizon as described in detail in [36]. 
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Given the energy functional E{v) we simply set 



^ '2G dv 



(6.16) 



(compare this with (6.8)). Then given this normalization, the dynamical horizon flux law 
becomes 



V2 



AT,,. 



iianb 
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(6.17) 



E{v2) - E{vi) = dv < j e 

Although this flux law is similar to (6.9), there are obvious differences. Specifically, the 
and Tabi'^n'' terms are no longer present — they have been absorbed into the definition 
of B, and consequently the scaling of the null vectors. In many ways, this is preferable as 
these terms appear to be associated to intrinsic features of the horizon and not to fluxes of 
matter or gravitational energy through the horizon. In particular, for a charged black hole 
e = Q^/R"^, while for the Kerr metric e is a function of J/M"^. 

Despite the similar forms and origins of the dynamical and slowly evolving flux laws, a 
closer analysis shows that they are actually quite different. We have managed to rewrite the 
dynamical horizon law in a manner which more closely resembles the slowly evolving horizon 
result. However, differences still remain. Further discussion of the dynamical horizon flux 
law and its relation to other formalisms may be found in [15, 17, 19, 23, 30]. 



VII. CONCLUSIONS AND OUTLOOK 

By concentrating on the foliating two-surfaces that are the constituent parts of future 
outer trapping horizons we have constructed a framework that encompasses isolated, slowly 
evolving, and dynamical horizons. The techniques used are similar in spirit to previous 
geometric analyses of horizons [18, 31], although our focus here has been on the deforma- 
tions of MOTS in the full four-dimensional spacetime, rather than a three-dimensional slice. 
Furthermore, we have seen that horizon evolutions and variations are both governed by the 
same underlying equations. Additionally, this set of equations is responsible for the various 
flux laws associated with these surfaces. Most of these results have been obtained previously 
by various authors [7-11, 15, 16, 18, 20, 28] using a variety of methods. The contribution of 
this paper is to rederive many of these properties of horizons using a common geometrical 
framework which highlights the connections between the various results. 

Thus we have seen that along with the freedom to refoliate an isolated FOTH is the con- 
comitant rigidity that prevents us from varying its three-dimensional structure. In contrast, 
the uniqueness of the foliation for a dynamical FOTH is the flip side of their well-known 
lack of rigidity against out-of-surface deformations. The freedom to vary dynamical FOTHs 
is also equivalent to fact that there is no unique method for evolving a given FOTS into a 
FOTH. Instead its evolution is many-fingered and any FOTS can be (locally) extended into 
many different FOTHs. 

This provides some insight into one of the most interesting open problems about these 
quasi-locally defined horizons. It is well-known that apparent horizons (in this case defined 
as the outermost marginally trapped surface on a spacelike slice S) can discontinuously 
jump during especially dramatic events such as black hole mergers. However, it is widely 
suspected (see for example [5, 20, 30]) that many of these discontinuities may result from the 
interweaving of the foliation with a continuous 6'(£) = surface that is sometimes a FOTH 
but at other times fails to satisfy one or both of 6'(„) < and 6n0[i) < 0. 
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Therefore, it becomes important to understand the circumstances under which a FOTH 
may end. We have seen that a FOTH may always be locally extended however, this certainly 
does not guarantee the existence of global extension. For example, in [20] it was seen that 
on horizons in Tolman-Bondi spacetimes, it is possible for to switch signs and become 
positive. At this point, the FOTH ends while a future inner trapping horizon begins. In 
these examples, there exists a continuous three-surface foliated by 6(^i) = cross- sect ions, 
though only part of that surface can be identified as a black hole horizon. A question 
deserving further investigation is then: what happens to FOTHs under finite extensions? 
More specifically, can a FOTH always be extended into an unending structure that is foliated 
by 6(^1) = surfaces? These extensions could fail if, for example, the repeated deformations 
broke the spacelike nature of the foliation two-surfaces or if there are circumstances under 
which the extension is open but bounded - that is H = U^S^ with v £ (^^1,^2) where Vi 
and/or V2 are finite. 

On a related note we motivated our horizon definitions by the notion that the interior 
of a black hole should consist of all points that lie on a trapped surface. However, we have 
seen that future outer trapping horizons are typically only the boundary of the trapped 
region associated to a particular slicing of the given spacetime. Another important question 
is then: what is the real boundary of the trapped region? It is widely believed [5, 16, 30] 
that this boundary corresponds to the event horizon in spacetimes where this structure is 
defined, however this has not been proved. A study of the finite extensions of FOTS where 
one tries to "push" them towards the event horizon would be one way of gaining insight into 
this problem. 

Returning to the results of this paper, a particularly interesting application of the de- 
formation rules comes in the detailed investigation of the near-equilibrium slowly evolving 
horizons, originally introduced in [26]. We characterized these as being almost-null in the 
sense that the equations governing their evolutions are almost the same as those for null 
surfaces in general and isolated horizons in particular. Slowly evolving horizons were shown 
to obey a dynamical version of the first law of black hole mechanics, just as in standard 
thermodynamics near-equilibrium systems obey the TdS-iorm of the first law. This result 
was first reported in [26] but the derivations appear here for the first time. We also saw 
that this result follows from the same "constraint" equations that are responsible for the 
corresponding Hawking-Hartle formula for event horizons. The fact that both of these re- 
sults depend on the horizons being almost non-expanding (as well as some of the examples 
considered in [27]) leads one to speculate that an eternally slowly evolving horizon might 
be indistinguishable from the corresponding event horizon to the same order of accuracy for 
which the other results hold. Further the rigidity of isolated horizons against deformations 
suggests that slowly evolving horizons might also be rigid up to this order. 

Our focus on the deformations of two-surface has been demonstrated to unify and illumi- 
nate diverse results in the study of quasi-local horizons. It seems likely that it will continue 
to be useful in studying future problems, including those outhned above. 
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APPENDIX A: DERIVING THE TWO-SURFACE EQUATIONS 



In this Appendix we catalogue some useful relations for two-surfaces embedded in four- 
space and use them to derive the equations of section II. 

1. Preliminaries 

It is often useful to decompose the four- dimensional Riemann tensor into its Weyl and 
Ricci tensor components. Thus we note that [1]: 

T^abcd = Cabcd + {ga[J^d]b + 5'fe[rf'^c]a) " g {ga[c9d]b) ^ , (Al) 

where Cabcd is the Weyl tensor and as usual square brackets indicate anti-symmetrization. 

Next, given a two-surface {S,qab,da) in four-space {M, gab,Va) the Gauss, Codazzi, and 
Ricci equations relate the curvature of spacetime to the intrinsic and extrinsic geometry 
of the two-surface. They may be derived fairly directly from a few facts. First, the push- 
forward of the inverse two-metric on S can be written as 

~ab ^ gab ^ ^a^b ^ ^b^a ^ ^^2) 

for any two null normals to 5* that satisfy £"77,^ = —1. Second, covariant derivatives of tensors 
defined intrinsically to S may be written in terms of the full four-dimensional derivative with 
the help of appropriate projections. Thus, for example, for a one-form iia G T*5', 

= Qa&'^di&f'^gfih) . (A3) 

Finally, the Riemann tensors on M and S respectively satisfy 

(VaVfc - V5Va)l^c = TlabcdW and (A4) 

{dadb - dbda)Wc = RabcdW'^ , (A5) 

for one-forms Wa G T*M and Wa G T*S. 

Then substituting (A5) into (A3) and doing some algebra, one can show that 

ra'll&d'J^ef,, = Rabcd + (€^^1^ + ^i^^^) ' (k^'J k^:J + k^^J kf,) , (A6) 

where k^^^ and k^^^ are the extrinsic curvatures defined in (2.2). This is the Gauss relation. 
Through similar calculations one can also derive the Codazzi relations: 

rjq^tn^fgh = {da - Ua)k^} " {d^ - Ub)k'^} (AT) 
faiqln^-^efgh = (da + CUa)ki:^ " (4 + CUb)kt^ 

where Ua is the normal bundle connection defined in equation (2.9). Alternatively, applying 
(2.3) and (Al) these may be expanded as 

{da - Ua)9^i) = 2(4 - Ub)a^^^ - ctaGbct - ^tCbcde^ I'^rf and (A8) 
{da + cDje(„) = 2(4 + c:;,)a(")^ - q\Gbjf + 2tCbcden''i''rf . (A9) 

Finally, by the same kinds of calculations used to derived the Gauss and Codazzi relations 
we can also derive the Ricci relation 

&n^C,,,^ = q:qttnfn,,,f = Qab - + ^i^j^^^ (AlO) 

where Qab is the curvature of the normal bundle defined in equation (2.10). 
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2. Deformation equations 

Next we consider the derivation of the deformation equations. As an initial step towards 
calculating Sx6(e) we find Sxh^^- To this end a few carefully chosen lines of algebra along 
with an application of equation (2.19) show that: 



(All) 



where the round index brackets indicate the usual symmetrization of the enclosed indices. 
Now equations (Al) and (A6) can be used to rewrite the first term, while the second can be 
shown to be 



Then, we find that 

^xkfh = -dadbB + 2iUi^adb)B + Kxk^ab 



(A12) 



(A13) 



+A [k^^lk'i^' - q'J%i^C,def - \qabGcd'i''^ 



-\-B ( -Rab 



1 



+B [ -^qlqt'R-cd + d(^aUJb) - tOaUJb 



With the help of (2.17) it is then straightforward to show that 



-d^B + 2Cj''d„B - B 



daCu'^ - R/2 + Gab^n' - e^e)Oin) 



-A [lla^f + G,,r£^+(l/2)^2j 



(A14) 



where HcjIP = uja^"" and ||cr(^)|p = cr^^V''^'*'^^. Similarly, we can take the trace-free part of the 
above to obtain 



ab 



-d{Jb}B + 2Cj{adb}B - Aqll'^qllfCdef 



(A15) 



+B 
+A 



^Ouwt} - ^0{n)(Tab + d{a^b} - ^{a^b} ' ^q{aqt}Gcd 



'mab 



1 

-( 

2 



1 

2' 



0", 



B 



"ac"b '"bc"a 



Here curly brackets around a pair of indices indicates their symmetric trace-free part (with 
respect to the two-metric). Thus, for example. 



^{adb}B = U[adb)B — -qaboj^dcB . 



(A16) 



Note that even though is trace-free, its variation will not usually inherit this property 
(this follows from the fact that 5x<f'^ 7^ 0). Thus, in the above expression, the first two lines 
are trace-free while the last line is not. 
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A direct expansion of q^C x{ncVb^'^) with applications of (2.19) and (Al) gives us the 
variation of the angular momentum one-form: 



Sx^a - daKx = -K> [d'A + Co'A] - k^^ [d' B - u'B] 

'I 



(A17) 



Finally, we can calculate the deformation of the two- curvature. Taking the standard varia- 
tion of the Ricci scalar (which is used, for example, in deriving the Einstein equations from 
the Einstein-Hilbert action [1]) and adapting it to two- dimensions, the variation 



5{m) = eiRab - l/2Rqab)6t + ed^ {r''t[dMbd - dd^q^c]) ■ 
Taking the 5s as 6xs and doing a little algebra this becomes 



6x{eR) = 2ed^d'{Aa^^ - Ba''^') -ed^AO^e^ - BO 



(A18) 



(A19) 



since the Rab = hRqab in two dimensions. Then, 



6xR = 2d''d\Aa^i - Ba'2') - d^AOi^,) - BO^r^)) - (A^(,) - 



(A20) 



APPENDIX B: MAXIMUM AND MINIMUM PRINCIPLES ON 5^ 

In this Appendix we briefly review the maximum principle for linear second order elliptic 
partial differential operators and then apply it to operators and functions defined on a surface 
which is diffeomorphic to S"^. 

As motivation we begin with a local maximum principle. Let U C IR** be an open set 
parameterized by coordinates where i G {1,2, .. .n}. Then any second order differential 
operator on the set of twice-differentiable functions over U takes the form 

Lm = a'^' . +/3^^ + 7/, (Bl) 

for some functions a*-', and 7 where i,j are indices and we assume the usual summation 
convention. If a*-' is positive definite then we say that L is elliptic. 

Now if / satisfies L[f] > everywhere on U for an elliptic L with 7 < 0, it is easy 
to show that / cannot have a non-negative local maximum in U. To see this recall from 
elementary calculus that if / has a local maximum at p, then df/dx'^ = for all i and the 
matrix d'^f/dx'^dx^ is negative definite. Equally elementary linear algebra tells us that for 
positive definite a*-' we must have 

d'^ f 

dx^dx^ ^ ^ 

at p. Then a term-by-term analysis of the right-hand side of (Bl) quickly shows that if 
/ > at p then < in contradiction to the original assumption. Therefore, if a local 
maximum exists, it must be negative. 

Note that as formulated this result doesn't cover cases where det{d'^ f /dx^dx^) = at 
fmax- However, the principle may be extended to all maxima and that is the content of the 
following theorem which is stated and proved in, for example, [40]. 
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Theorem 1 (Hopf's Maximum Principle) Consider a second order differential opera- 
tor 

on a connected open set U G H"" over which 7 < 0. Assume that the functions /3* and 7 are 
locally bounded and that in a neighbourhood of any point of U there are constants Mi and 
M2 such that 

Mi(e-o<«^^a, <M2(e-o, (B4) 

for all G ]R", where ■ is the usual Euclidean dot product for H" and we understand 
summation over repeated indices. Finally let f G C'^{U) and assume that it satisfies 

L[f]>0, (B5) 

everywhere in U. Then f cannot have a non-negative maximum on U, unless f is everywhere 
constant. 

We can then apply this theorem to show the following: 

Corollary 1 (Maximum principle on a two-sphere) Let S be a two-manifold that is 
topologically and has spacelike two-metric Qab and compatible covariant derivative da- 
Further, let f G C^(S') be a scalar field that satisfies 

L[f] > , (B6) 

everywhere on S for a differential operator of the form 

L[f] = r'dad,f + /3^dJ + ^f, (B7) 

where j3°' G TS and 7 < 0. Then f is either constant or everywhere negative. 

To see this first note that 5*^ is compact and so / must have (and achieve) an absolute 
maximum. However, if we consider any finite set of charts that cover S, Theorem 1 applies 
to the coordinate realization of L[f] on each of those charts. Thus working chart-by-chart 
and piecing the results together, either / is constant over 5* or the absolute maximum 
fmax < 0. That is, / is either everywhere constant or everywhere negative. □ 
Similarly with a simple substitution / ^ — / we have a minimum principle: 

Corollary 2 (Minimum principle on a two-sphere) Let S be a two-manifold that is 
topologically and has spacelike two-metric (jab and compatible covariant derivative da. 
Further, let f G C^{S) be a scalar field that satisfies 

L[f] > , (B8) 

everywhere on S for a differential operator of the form 

L[f] = -tdadbf + ^''daf + 7/ , (B9) 

where G TS and 7 > 0. Then f is either constant or everywhere positive. 

These results are used repeatedly in section III. 
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FIG. 4: The FOTS defining quantities for a Kerr liorizon. That < everywhere on all Kerr 
horizons is shown by a) which graphs it for the standard scaling of the null vectors. However b) 
and c) demonstrate that the 5n^(£) < condition is not so trivial. For the standard (C2) scaling 
of the null vectors SnO^i) can be positive and we must rescale as in equation (C4) to show that the 
condition is always satisfied. 

APPENDIX C: CLASSIFICATION OF KERR HORIZONS 

To justify some of the assumptions made in defining FOTHs and slowly evolving horizons, 
we consider the values that some of the key geometric quantities take in the Kerr spacetime'^. 
This should certainly contain both a FOTH and a slowly evolving horizon (with e = 0) so we 
must verify that the various assumption hold on the Kerr horizon. Further, for spacetimes 
that are perturbations of Kerr the values of these quantities should be similar and so these 
calculations gives us some intuition about these spacetimes and their classification as well. 

In generalized ingoing Eddington-Finkelstein coordinates, the Kerr metric takes the form 



1 



2mr 



dv + 2dvdr 



Aamr sin 6 



dvdcf) + T.de'^ 



(CI) 



sin^ 9 {{r'^ + a 



2\2 



sin^ 



2mr + a?)) 



where E = + cos^ 9 and the event horizon is at r = m + y'lm? — o?. 

It is common to scale the null vectors so that ^°^ is proportional to the global "time 
translation" Killing vector field on H and the corresponding flow evolves the v = constant 
two-surfaces into each other. Identifying this scaling with a sub/superscript-"i^'" we have 



-K 



n 



K 



1,0,0, 



+ 



and 



(C2) 



sin^ 9 
2 + a^) 



^ Similar results hold for the Kerr-Ncwmann family, but for ease of presentation we restrict our attention 
to pure Kerr. 
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FIG. 5: The quantities that determine SnO^e) (which also appear in the first condition of the 
definition of a slowly expanding horizon) as calculated for the standard scaling of the null vectors 
(C2). 

where the coordinate ordering is {v,r,6,(f)}. This is also the scahng for which the surface 
gravity is constant and takes the value 

In Figures 4a and 4b we plot 6'(„) and SnO^i) respectively for this scaling of the null 
vectors. It is immediately apparent that while 6'(„) < 0, horizons with sufficiently large 
angular momentum have SnO{£) > over some region. However, this does not imply that for 
these horizons the v = constant surfaces fail to be FOTHs. Instead it nicely demonstrates 
that different scalings of the null vectors generate different deformations of the two-surfaces 
and that not all of the resulting two-surfaces are fully trapped. If we rescale the null vectors 
to become 

^TH = ^^K and = Snf , (C4) 

then SnO{e) calculated with respect to these vectors is everywhere negative, as shown in 
Figure 4c. Thus the v = constant slices are FOTS, but the standard scaling of the Kerr null 
normals does not show this. 

Computationally we can understand what is happening for these two different scalings 
by considering the individual components that went into the calculations. Equation (3.1) 
tells us that in the absence of matter fields 

SJ^^^ = -R/2 + Ipf - dau" . (C5) 

Figure 5 graphs these quantities for the standard scaling (C2) of the null vectors. The Ricci 
scalar is, of course, scaling invariant and Figure 5a shows that for horizons with large angular 
momentum it can be negative. Thus, the scaling dependent daUi"" term must be positive and 
large enough to compensate for this. However, a quick examination of the shapes of Figures 
5b and 5c shows that for the standard rescaling, this term is also negative. 

In fact these figures are also sufficient to show why the rescaling (C4) resolves this problem 
— although in a slightly roundabout way. To see this we first note that for any (topologically) 
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a) 



a/ liH 



b) 



FIG. 6: The quantities appearing in the definition of the slowly evolving parameter e as calculated 
for the standard scaling of the null vectors (C2). 

spherical surface embedded in a spacetime, there is always a scaling of the null vectors so 
that daOj"" = (as discussed in [12, 41] this ultimately follows from the Hodge decomposition 
theorem). For the two-surface under consideration this scaling is 

C = -^tK and n° = v^nf . (C6) 

Taking this as a reference, equation (2.9) tells us that if = and Ua = {l/f)n° then 

||cDf = llcDof + 2cuXln/ + ||rfln/f and (C7) 
daC^" = dHnf . (C8) 

On an axisymmetric horizon, the second term on the right-hand side of (C7) will vanish. 
Meanwhile, ||u;|p is invariant under a rescaling / — ^ 1//, while daO"' —daCb"-. Just such a 
rescaling is made in the change from the Killing normalized null vectors (C2) to the FOTH 
normalized ones (C4), where in this case / = 1/a/S. Thus, it follows that 

= -R/2 + \\uKf + 4cD^ , (C9) 

and in this case the daio"' term is preciesly what is required to guarantee that 5n0{^ti > over 
the whole horizon. 

Figure 5 also demonstrates that R, and daCo"" are each of order l/Rfj as required 

for a slowly expanding horizon. Combining this with the fact that the expansion parameter 
C = for the v = constant two- surfaces, we confirm that the event horizon of a Kerr black 
hole is not only a FOTH but is also a slowly evolving horizon. 

Finally, Figure 6 graphs the quantities that appear in (5.6) in the definition of a slowly 
expanding horizon. Although in this case their exact form doesn't matter (since C = e = 0), 
it is useful to note that they too are of order l/Rfj- Thus, for perturbed Kerr solutions (5.6) 
will be satisfied and the horizon will be slowly expanding, as long as C is sufficiently small. 
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